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The cold dark matter (DM) paradigm describes the large-scale structure of the universe remark- 
ably well. However, there exists some tension with the observed abundances and internal density 
structures of both field dwarf galaxies and galactic satellites. Here, we demonstrate that a sim- 
ple class of DM models may offer a viable solution to all of these problems simultaneously. Their 
key phenomenological properties are velocity-dependent self-interactions mediated by a light vector 
messenger and thermal production with much later kinetic decoupling than in the standard case. 



PACS numbers: 95.35. +d, 95.85.Pw, 98.65.-r, 98.70.Sa 

Introduction. — Recent advances of cosmological pre- 
cision tests further consolidate the 'cosmological concor- 
dance model', indicating that 4.5% of the mass in the uni- 
verse is in baryons, 22.6% is non-baryonic cold dark mat- 
ter (CDM), and the rest is Einstein's cosmological con- 
stant A (or behaves like it) pQ. The leading CDM candi- 
dates are weakly interacting massive particles (WIMPs) 
which are thermally produced in the early universe [2J. 
While their chemical decoupling from the heat bath sets 
the observed DM relic density today, their kinetic decou- 
pling induces a small-scale cutoff in the primordial power 
spectrum of density perturbations [3J- For neutralino 
DM, e.g., this cutoff corresponds to a smallest protohalo 
mass of M cut /M© - Kr u -l(r 3 @J, but it could be as 
large as M cut > 10 M if DM couples to new light scalars 
0. After kinetic decoupling, standard WIMP CDM be- 
haves like a collision-less gas. Baryons, on the other hand, 
can radiate away excess energy and sink to the centers of 
CDM halos where they form stars and galaxies. In this 
picture, structure formation proceeds hierarchically with 
galaxies to form at sites of constructive interference of 
small-scale waves in the primordial density fluctuations. 

Despite the great success of ACDM cosmology, detailed 
observations of nearby small galaxies pose a number of 
puzzles to this paradigm. Here, we isolate three distinct 
classes of problems. (1) The observed galaxy luminos- 
ity and Hi-mass functions show much shallower faint-end 
slopes than predicted by ACDM models [5] ; this is locally 
known as the 'missing satellites problem' of the Milky 
Way (MW) , which should contain many more dwarf-sized 
subhalos than observed [7J. (2) Simulations predict an in- 
ner DM cusp for the density structure of galaxies, seem- 
ingly at odds with the cored profiles found in observed low 
surface brightness galaxies and dwarf satellites [5]. (3) 
Recently, it was realized that the most massive subhalos 
in ACDM simulations of MW-size halos have an internal 
density structure that is too concentrated in comparison 



to the observed brightest MW satellites: the simulated 
circular velocity profiles increase more steeply and attain 
their maximum circular velocity at smaller radii than any 
of the observed ones. On the other hand, those simu- 
lated subhalos should be 'too big to fail' in forming stars 
according to our understanding of galaxy formation (be- 
ing more massive than the UV-photosuppression scale at 
all redshifts, after formation, for conceivable reionization 
histories). Thus, it is extremely puzzling why there is no 
observed analogue to those objects [9]. 

Astrophysical solutions to (1) invoke suppressing the 
formation of galaxies within existing dwarf halos or sup- 
pressing the star formation in dwarf galaxies. Galaxy 
formation can be held back by increasing the gas en- 
tropy before collapse, e.g. via photoionization |10) . blazar 
heating [TT] or AGN feedback in the radio-quiet mode 
|12) . A photoionization-induced lack of Hi [T3] or in- 
trinsically low metallicities |14j may further suppress the 
cooling efficiency of collapsing baryons. Numerical sim- 
ulations with a photoionizing background, however, can- 
not suppress dwarf galaxy formation at the level implied 
by observations [15j . In principle, gas may also be re- 
moved from dwarfs via photo-evaporation [16 and feed- 
back from supernovae [17j . Any such feedback, however, 
implies remnant stellar populations and Hi masses in con- 
flict with most recent observational constraints |18j . 

The 'cusp-core' problem (2) may be addressed by large 
velocity anisotropies or reduced central DM densities. 
There is a degeneracy between cored isotropic and cuspy 
anisotropic velocity distributions and the stellar line-of- 
sight velocity data is still too sparse to dynamically re- 
solve (2) [19]. Reducing central DM densities was pro- 
posed as a result of efficient baryonic feedback processes 
|20) . however in contradiction to cuspy dwarf profiles in 
other simulations with feedback |21j . 

The 'too big to fail' problem (3) might be solved by 
either an increased stochasticity of galaxy formation on 
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these scales or a total MW mass < 8 x 10 11 M Q [2"2"] . 
Abundance matching of stellar and halo mass, which 
agree with stacking analyses of gravitational lensing sig- 
nals and satellite dynamics of SDSS galaxies, make the 
required large degree of stochasticity implausible [2"5] . 
For a 10 12 M MW, on the other hand, the chance to 
host two satellites as massive as the Magellanic Clouds is 
less than 10% [21] and even lower for smaller MW masses 
(from satellites studies of MW-type SDSS systems [25] 
and MW and Andromeda orbit timing arguments |26j). 

The next logical possibility that could lead to a sup- 
pression of small-scale power is a modification of the 
CDM paradigm itself. The most often discussed options 
are interacting DM (IDM) ;27] and warm DM (WDM) 
[2"8] , though it should be noted that there exist interesting 
alternatives such as DM from late decays [23], DM with 
large annihilation rates [3U] , extremely light DM particles 
forming a condensate [31] . or inflationary models with 
broken-scale invariance [32] . As was soon realized, how- 
ever, IDM with a constant cross section produces spheri- 
cal cores in conflict with observed ellipticities in clusters 
[3"3"] and the survivability of satellite halos [31] • While 
WDM is unlikely to account for some of the large ~1 kpc 
cores claimed in dwarfs [35] , and severely constrained by 
Lyman-a observations 36435]. it may be able to par- 
tially resolve the 'too big too fail' problem by allowing 
these subhalos to initially form with lower concentrations 
[ID] . Alternatively, DM self- interactions mediated by a 
Yukawa potential, with the resulting characteristic veloc- 
ity dependence of the transfer cross section [UJ[12], avoid 
constraints on scales of MW-type galaxies and beyond 
[4"3"] and produce ~lkpc cores that match the observed 
velocity profiles of massive MW satellites [2] (see also 
Ref. [31]). 

Most astrophysical and DM solutions have shortcom- 
ings, or can explain at most two of the three problems, 
which makes them less attractive on the basis of Oc- 
cam's razor. Here, we demonstrate that there is a class of 
IDM models that simultaneously can account for all three 
problems. Encouraged by the results of Refs. [UJUl], in 
particular, we will focus on models with a Yukawa-like 
interaction between the DM particles that is mediated 
by a light messenger (see Fig. [I]). As we will show, the 
kinetic decoupling of DM in these models can happen 
sufficiently late to suppress the power spectrum at scales 
as large as that of dwarf galaxies, M cut > 10 9 M Q , while 
at the same time the velocity-dependent self-interaction 
of DM produces cored density profiles in dwarfs 1 . 

Model setup. — In models with new light exchange 
particles <j), kinetic decoupling can happen much later 



Late kinetic decoupling (T^j ~ 0.1 keV) was advocated as poten- 
tial solution to the missing satellites problem before |45| - though 
this analysis considerably under-estimated for WIMPs |46| . 



than in standard WIMP scenarios, in particular for small 
masses m<p [5j. For scalar exchange particles, however, 
the amplitude for DM scattering with leptons scales like 
~ m x mi/m^ ) , implying that scattering with neutrinos is 
generally negligible. While a coupling of <fi to charged 
leptons also leads to a loop-suppressed effective coupling 
to photons, C D g^cfiF^ F^, the resulting scatter- 
ing amplitude does not contribute in the relevant limit 
of small momentum transfer. Kinetic decoupling there- 
fore never occurs at <C 0.1 MeV, at which point the 
number density of electrons starts to become strongly 
Boltzmann-suppressed and there are no lighter (and thus 
more abundant) particle species left that could keep up 
kinetic equilibrium instead. 

Let us consider instead the situation where DM con- 
sists of heavy Dirac fermions \ which only couple to a 
light vector boson V. Due to our interest in late kinetic 
decoupling, we will require V to also couple to neutrinos: 

Ant => ~g x xfx ~ g v vYv . (1) 

Note that we take a phenomenological approach here and 
only state couplings that explicitly enter our analysis. In 
particular, V does not have to be a gauge boson, which 
leaves couplings to other SM particles unspecified (see 
e.g. Ref. [13 for a recent model- independent analysis). 
DM is then thermally produced in the early universe via 
XX VV. Assuming g v is small, but large enough to 
thermalize V at early times, the relic density is given by 




This expression receives 0(1) corrections due to the Som- 
merfeld effect [48], i.e. a multiple exchange of V as shown 
in Fig. [T] which we fully take into account in our anal- 
ysis. The kinetic decoupling temperature, on the other 
hand, will be set by \~ v scattering. The corresponding 
amplitude at small momentum transfer reads 

all spins ^ 

In the following, we will consider g v as an essentially 
free parameter while g x is fixed by the requirement to 
obtain the correct relic density (see e.g. Ref. [49 for a 
list of possible natural explanations for g v <C g x ). 

DM self- scattering. — The light vector messenger in- 
duces a long-range attractive Yukawa potential be- 
tween the DM particles, cf. Fig. [T] Concerning elas- 
tic DM self-scattering, this is completely analogous to 
screened Coulomb scattering in a plasma for which sim- 
ple parametrizations of the transfer cross section ctt(w) in 
terms of m x , my , g x and the relative velocity v of the DM 
particles exist [HI [SDJ ■ Using these parametrizations, it 
was shown that the type of DM model introduced above 
produces cores rather than cusps ,43] and may solve the 
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not enough flattening 
of cuspy profiles 



Figure 1: Interaction processes that set the DM relic density 
and may lead to observable neutrino annihilation products 
today (left), change the inner velocity and density profile of 
dwarf halos (middle) and induce a comparatively large cutoff 
in the spectrum of primordial density perturbations (right). 



'too big too fail problem' [33], without being in conflict 
with the strong constraints for models with constant o~t- 
We also note that <7t drops with larger v such that for 
galaxy clusters only the very central density profile at 
r < 0(1 — 10)kpc will be smoothed out, matching ob- 
servational evidence (from improved lensing and stellar 
kinematic data [ST]) for a density cusp in A383 that is 
slightly shallower than expected for standard CDM. 

For our discussion, the astrophysically important 
quantities are the velocity f^ax = 9 x m v /(2Tr 2 m x ) at 
which utv becomes maximal and <r™ ax = <7T(^max) = 
22.7 my 2 . In particular, u max should not be too differ- 
ent from the typical velocity dispersion o~ v ~ O(10) km/s 
encountered in dwarf galaxies if one wants to make any 
contact to potential problems with standard structure 
formation at these scales. On the other hand, the value 
of <7™ ax is constrained by various astrophysical measure- 
ments, see Ref. [44; for a compilation of current bounds. 

Fixing g x by the relic density requirement, there is a 
one-to-one correspondence between the particle physics 
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Figure 2: The white area corresponds to DM and mediator 
masses that may solve the 'cusp vs. core' problem. The crosses 
indicate two benchmark models for which detailed simulations 
[44] have found a solution to the 'too big to fail' problem. 
Dashed and solid lines show contours of the astrophysical rel- 
evant quantities u^ ax and tWx- See text for further details. 
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input (m x ,my) and the astrophysically relevant param- 
eters (wmax, cr™ ax ). As demonstrated in Fig. [2j a so- 
lution to the aforementioned small-scale problems (2) 
and (3) may then indeed be possible for DM masses 
of m x > 600 GeV and a mediator mass in the (sub-) oTta^aU^ 



Figure 3: This plane shows the mediator mass my vs. the 
coupling strength g v . Large values of g v and small values of 
my lead to late kinetic decoupling and thus a large mass M cut 



MeV range. We also display the strongest astrophysi- 
cal bounds on large DM self-interaction rates [43]. For 
m x ~ ^ TeV, they arise from collisions with particles from 
the dwarf parent halo, while at larger m x an imminent 
gravothermal catastrophe is more constraining. 

The small-scale cutoff. — For small kinetic decou- 
pling temperatures Tkd, acoustic oscillations [52] are 
more efficient than free streaming effects to suppress the 
power spectrum [JJ I53j . The resulting exponential cutoff 
can be translated into a smallest protohalo mass of 



by Ly-a data while M cut > 1O 9 M0 may solve the small-scale 
abundance problems of ACDM cosmology. 



those expressions to allow for T v ^ T, we find 
0.062 keV (T 



NJ (g x 9u) 1 



kd 



4tt 

3 m 



= 1.7 x 10' 



\keV 



T=T k , 



Mr. 



(4) 



where H is the Hubble rate and we assumed late kinetic 
decoupling such that the effective number of relativistic 
degrees of freedom g e f[ = 3.37. For scattering with rela- 
tivistic neutrinos, c.f. Eq. ([3]), the analytic treatment of 
kinetic decoupling given in Ref. |54j is valid. Extending 



where N v is the number of neutrino species coupling to 
V . Combining this with Eq. ^ we therefore expect that 
T^d, and thus M cut , is essentially independent of g x and 
m x . 

Using for definiteness N u — 3 and T v = (4/11) s T 7 , we 
show in Fig. [3] contours of constant M cut in the (g v ,mv) 
plane. We find that the result of the full numerical 
calcul atio n [31 [S] is indeed extremely well described by 
Eqs. ( 4p I for g v > 10~ 7 (assuming m x ~ 1 TeV and 
my ~ 1 MeV; this value is even lower for larger m x and 
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smaller my). For g v < 10 -7 , DM scattering with the 
non-relativistic mediator particles V starts to dominate 
over scattering with neutrinos, and M cut eventually be- 
comes independent of g v . We checked that the new era 
of DM annihilation generally expected in models with 
Sommerfeld-enhanced annihilation rates (see Ref. [5] for 
a consistent treatment) has only negligible impact on our 
results for the late decoupling times we focus on here. 

Lyman-a forest bounds. — Conventionally, a possi- 
ble cutoff in the power spectrum is often expressed in 
terms of the mass m w of a WDM thermal relic. In this 
case, it is set by free streaming of the WDM particles 
and the comoving free-streaming length Rf is given by 

Rf = 0.1 (O m ^ 2 /0.13) 1/3 (m w /keVr 4/3 Mpc 55J. For 
a characteristic wavenumber kf = 0.46/i?f , the linear per- 
turbation amplitude is suppressed by a factor of 2 and 
the characteristic filtering mass can be defined as [55] 

„, a f^) S = 5 , xl ».» er V., (6) 

where Af = 2-njkf ~ 13.6i?f. This choice of Mf is justi- 
fied by numerical experiments [56j that find the resulting 
halo statistics for an initial density distribution with a 
sharp cut in the power spectrum at k c — 27r/Af to be 
very similar to the statistics of an initial density field 
smoothed with a top-hat window of radius Af/2. Cosmo- 
logical WDM simulations show a deviation of the mass 
function from the CDM case on scales given by Eq. © 

0[57]. 

Combining data of the Lyman-a forest, the cosmic mi- 
crowave background (CMB) and galaxy clustering allows 
to constrain the cutoff scale in the power spectrum; in 
terms of the mass of a thermal WDM candidate, a 2a- 
bound of m w > 2 keV has been claimed [35J ISZ] ■ This 
weakens to m w > 0.9 keV when rejecting less reliable 
data at z > 3.2 [35] due to systematic errors [35]. Revisit- 
ing Lyman-a data yielded m w > 1.7keV which, however, 
is subject to systematic uncertainties at the ~30% level 
[39], especially considering that blazar heating was not 
accounted for in deriving cosmological constraints [58j . 

Lyman-a data thus firmly exclude m w < 1 keV or 
Mf > 5.1 x 10 10 M (corresponding to a maximal cir- 
cular velocity u max ~ 70 km s _1 ). For m w ~ 1 — 2 keV, 
WDM models are able to alleviate the 'missing satel- 
lite problem' (somewhat depending on feedback recipes) 
[59"] . bring the faint-end of the galaxy luminosity func- 
tion into agreement with data [60) . and match [5] the 
Hi velocity function measured in the ALFALFA survey 
2 . For m w > 3 keV, the corresponding mass cutoff 
Mf < 6 x 10 s M Q is too small to have any impact on the 



2 The small-end of simulated velocity functions for a thermal relic 
WDM particle with m w < 3 keV are noticeably shallower than 
the analytical estimates of the Sheth and Tormen [61 1 - 



faint-end of the galaxy luminosity function. We include 
these bounds in Fig. [3] to demonstrate that our model 
can successfully address also the abundance problem (1). 

Discussion. — In a phenomenological approach to 
identify the key properties of DM models that can ad- 
dress all three ACDM small-scale problems simultane- 
ously, we found it sufficient to simply postulate the exis- 
tence of a light vector messenger V that couples to both 
DM and neutrinos as in Eq. ([!]). If V does not couple 
to quarks or other leptons, the coupling g v is essentially 
unconstrained [47]. While beyond the scope of this let- 
ter, however, we stress that it would be very worthwhile 
to study possible concrete realizations of our setup. 

The greatest challenge for such a model building might 
be to prevent, even at the one-loop level, a kinetic mixing 
between V and photons, which is severely constrained for 
m v ;$ MeV [53]. On the other hand, limits on tree-level 
couplings of V to charged leptons (e.g. from contributions 
to the anomalous magnetic moment of the muon, beam 
dump experiments or low-|<7| 2 v-e scattering [62] [63] ) 
seem less severe and could at least partially be evaded 
by generation-specific couplings. Another option could 
be a new U(l) coupling to DM and sterile neutrinos v s 
|64) . As long as the v s have been in equilibrium in the 
very early universe and are relativistic at Tkd, this would 
not change the phenomenology of our model; Eq. ^ , in 
particular, would still apply. It is therefore quite inter- 
esting that CMB observations seem to favor additional 
relativistic degrees of freedom, which corresponds to the 
presence of one light sterile neutrino species [55] . 

Finally, we note that for typical galactic velocities 
v ~ 10 -3 , the type of DM candidate we propose here 
annihilates with a Sommerfeld-enhanced rate of (av) ~ 
3 x 10~ 24 (m x /TeV) _2 cm 3 s _1 into a VV pair which then 
decays exclusively into neutrinos (if my < 2m e ). Such a 
large annihilation rate will be in reach of future IceCube 
observations of the galactic center [55] ; for m x < 1 TeV, 
in fact, a strong Sommerfeld-induced substructure en- 
hancement of the signal [57. may already be constrained. 

Conclusions. — We have introduced a class of DM 
models with the unique property of addressing all 
three ACDM small-scale problems simultaneously, which 
should make them very attractive alternatives to be stud- 
ied. From a model-building point of view, the only in- 
gredient that is needed is a (sub-)MeV vector messenger 
particle that weakly couples to neutrinos and even more 
weakly to other standard model particles. While collider 
and direct searches for DM will be extremely challenging 
in this scenario, a TeV neutrino signal from the galactic 
center could turn out to be a smoking gun signature. 
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